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Abstract 

We analytically compute asymptotic expansions of a 1-dimensional sub- 
manifold of stable and unstable manifolds in a 4-dimensional symplectic 
mapping by using the method called asymptotic expansions beyond all or- 
ders. This method enables us to capture exponentially small splitting of 
separatrices and also to obtain explicit functional approximations of the 
sub-manifolds. In addition, we show the condition with which homoclinic 
structure caused by crossing between the stable and unstable sub-manifolds 
is regarded as a direct product of 2-dimensional mappings. 
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1 Introduction 



Existence of transversal homoclinic (or heteroclinic) points implies non-integrability 
of Hamiltonian systems and it is important to understand homoclinic structure 
for studying phase space's structure of Hamiltonian systems. Melnikov's method 
(or Melnikov's integral) is a powerful tool to search for transversal homoclinic 
points This method is based on the regular perturbation method with respect 
to perturbation strength e, hence it is difficult to apply this method to singular 
perturbed problems straightforwardly. 

However in Hamiltonian systems with rapidly forces, and 2-dimensional sym- 
plectic mappings which are perturbations from the identical mappings, and so on, 
it is known that the homoclinic structure contains singular terms, for example, 
exp {—P/e"') ,a,P > 0-^, and concequently had not been captured with the 
regular perturbation method with respect to e. 

The difficulty however have been overcome by using the method called asymp- 
totic expansions beyond all orders. This method was first proposed and applied 
to the standard mapping by Lazutkin and co-workers[^. They calculated inter- 
secting angles between stable and unstable manifolds asymptotically. And the 
method was improved by using Borel transformation and Stokes phenomenon [^, 
1^. The same approaches were developed independently by Kruskal and Segur[Q 
in the model of crystal growth and by Pomeau et al. in the K-dV solitons under 
a 5-th order singular perturbation. Afterwards Tovbis et al. applied this method 
to Henon map [p!0|, p!T[] and constructed explicit functional approximations of the 
stable and unstable manifolds. Nakamura and Hamada analyzed a 2-dimensional 
symplectic mapping with a double well potential in the same way and calculated 
intersecting angles of stable and unstable manifolds [|T^ . 

The applications mentioned above are all to 2-dimensional symplectic map- 
pings, which correspond to Poincare mappings of Hamiltonian systems with two 
degrees of freedom. On the other hand, 4-dimensional symplectic mappings corre- 
sponds to Poincare mappings of systems with three degrees of freedom. It is very 
impotant to extend the analytical method for 2-dimensional symplectic mappings 
to 4- or more-dimensional symplectic mappings, because phase space's structure 
of Hamiltonian systems with three or more degrees of freedom quite differs from 
that of systems with two degrees of freedom, e.g., Arnold diffusion . 

An application to 4-dimensional symplectic mappings is reported by Gelfreich 
and Sharomov||13||. They treat a coupled 4-dimensional standard map and prove 
the splitting angle between 1-dimensional stable and unstable sub-manifolds is 
exponentially small and the splitting angle orthogonal to the sub-manifolds is 
not exponentially small. Their model however contains the particular coupling 
potential, i.e., two-body central force coupling J{qi — ^2) (compare the mapping 
® or (i). 

We would like to take more general 4-dimensional mappings into consideration 
and construct functional approximations of 1-dimensional sub-manifolds of stable 
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and unstable manifolds. In this paper we shall show that it is possible to apply 
this method to more general coupling potential terms. 

In this paper we consider a 4-dimensional symplectic map 

dJ 

Pj,n+i = Pj,n - e (2g|_„ - qj^n) - ^^^^'^'q^ ^ j = 2, (1) 

'lj,n+l = (lj,n + ^Pj,n+1 

where < e ^ 1, k G R, 7 G N, and J = qija^2%.^ ai, ^2 G N is a coupling term. 
The map (|l|) is an exact symplectic map, hence there exists a generating function 
W = I + tJ [|15l, where 

/ = Pj,n+iqj,n, (2) 

J = E{^fe-,"+i) + ng,,n)} + eWfe,„), (3) 

Tip) = \v\ (4) 
^(^) = (5) 

And the map (P is defined as gj>+i = dW jdpj^n+i and pj^n = dW/dqj^n- For 
e = the map is obviously identical, hence the map (|l|) can be regarded as 
a perturbation from the 4-dimensional identical map. Note that the map (|l|) 
becomes chaotic when e 7^ 0, in general. 

The symplectic map ([^) has a fully hyperbolic fixed point, which is a direct 
product of two of 2-dimensional hyperbolic fixed points, at the origin (0, 0, 0, 0) 
(see the equation (^). Hence there exist 2-dimensional stable and unstable man- 
ifolds near the origin. 

In the present work we take ai = 02 = 2, 7 = 2, i.e., J = qfnQ2n simplify 
the following computation and notation. Thus the symplectic map (|l|) is rewritten 
as 

3 ^ \ 0^3,^^ ^2 



\ Pj,n+l - Pj,n e(2g^> qj,n) 2e Kqj,nqk,n ^ j ^ ^ ^ 2} , j k. (6) 
I Qj,n+1 = qj,n + ^Pj,n+1 

To construct the unstable manifold we transform the map into the second 
order difference equations and change the independent variable n into t with 
meanings of time. We take e as the difference parameter, that is, qj,n±i = qj{t ± 
e), J = 1, 2 if g,- „ = q^{t) Thus 

A%it) = qj{t)-2q^{tf~2e^Kq^{t)qk{tf (7) 

where Alq,{t) = {q,{t + e) - 2g,(t) + - e)}/e\ j = 1, 2. 



2 



One can transform the equations (|^ into singular-perturbed ODEs|T0|-[T^, 

^ = g.(t)-2g,<tf-2e^.g,<t)g.(t)^-2£^^^ (8) 

One of the purposes of this paper is to construct explicit functional approxima- 
sions of 1-dimensional stable and unstable sub-manifolds which satisfy qi {t, e) = 

An outline of constructing functional approximations of the sub-manifold, 
which is also an outline of this paper, is as follows: We construct solutions to the 
equations (H) with the regular perturbation method (section Q). The solutions 
however break down around t ~ because of existence of singular points. Hence 
we blow up complex t-plane near one of the singular points and analyze with 
singular perturbation methods (section H). One can rewrite the equations (|^) 
near the singular points (section ^ as 

A2$j = -2<l>i^ + e2$j - 2Ke2$j$2, (9) 

where A'^^j{z) = ^j{z + 1) - 2^j{z) + ^j{z - 1), and j,k e {l,2},j 7^ k. 
We analyze the equations (|^) by using Borel transformation, which enables us to 
capture Stokes Phenomenon (section ^. We match the solutions of the equations 
(^ to the solutions of the equations (^ and successfully obtain the analytical 
representation of a 1-dimensional sub-manifold of an unstable manifold in the 4- 
dimensional symplectic mapping (section Finally we will conclude the result 
(section |^ and give a summary of this paper in section ^. 



2 The outer equation 

In this section we try to construct a 1-dimensional sub-manifold of an unstable 
manifold of the equations (^) by using the regular perturbation method, i.e., 
Taylor expansion with respect to e. Although the solution which satisfies the 
boundary conditions is obtained, it does not converge. We must therefore treat 
the equations (^) with the singular perturbation method in the next and following 
sections. Outline of construction of the manifold, which corresponds to contents 
from section ^ to section |^, goes along the way by Tovbis et al.[|l^, From 
now on, the subscripts j and k will be used as the meaning of j, k G {1, 2}, j 7^ k 
without further comment. 



First using Taylor expansions as 

00 

q,it±e) = Y: 



1=0 



n M ' 



one can rewrite the equations (0) as 



^ = g.(t)-2g,(tf-2e^.g,(t)g,(t)^-2£^^. (10) 
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Equations ( p!0[ ) and solutions to them are called the outer equations and the 
outer solutions, respectively. One can construct asymptotic representations of an 
unstable manifold as t — > —oo by solving the ODEs (|10|). Boundary conditions 
with which we solve the equations (0) are 



hHi^gj,„(t,e) = 
Qj^uit, e) > 



(11) 



where the subscript u stands for an unstable manifold. The conditions (pJ]) 
mean that the unstable manifold is asymptotic to the hyperbolic fixed point with 
positive value. 

We construct a formal solution qj^it, e) of the outer equations ([ToD in the 
power series of as 

oo 

koit,e) = T.1Jfii{ty'^ (12) 

1=0 



which we call the outer expansions. Substituting the equations ([l^ ) into the 
equations ([T0|), equations are successively obtained for each power of e^; 

O(,0) . ^ g,.oo(t)-2g,,oo(t)^ , (13) 

0{e'') : = {1 - QqMtfhMt) + F,,oi (^ > 1), (14) 

where Fj^i are polynomials of gm,oo(^), (lm,oi(t), ■ ■ qm,o,i-i(t),m = 1, 2 and their 
derivatives. The boundary conditions (0) are rewritten as 



lim qj.oiit) = ] 
qjfiiit) > J 

for / = 0, 1, ■ ■ •. 

From the equations ( |T3D and ([T5|) we obtain unperturbed solutions as 



(15) 



Qj,oo{t) = sech(t — t 



where tjS are integral constants. 

In this paper we restrict our attention to a 1-dimensional sub-manifold which 
satisfies gi^u(t, e) = q2,u{'t,^), because it seems that construction of the whole 
2-dimensional manifold is difficult. Hence one must put ti = ^2- Furthermore 
the system is autonomous; accordingly we can put ti = ^2 = without loss of 
generality. Thus 

qjfioi^) = sech(t). (16) 



4 



The unperturbed solutions (|I^) possess even symmetry. Hence let us suppose 
that solutions of the outer equations (|10]) also possess the even symmetry, i.e., 

qj,u{-t,t) = gj>(t,e), j = 1,2. (17) 

Because the stable sub-manifold qj^sit,^) is obtained as qj^sit,^) = Qj,u{—t,^), if 
one can obtain consistent solutions of the equations (|T0|) under the conditions 
([TTD, the stable and unstable sub- manifolds coincide perfectly, i.e., homoclinic 
bifurcation does not occur. Expanding the conditions (0) in the power of one 
obtains 



qjfii{-t) = (imit)^ / = o,i,---. (18) 

Next, by using the equations (|1^) we shall construct solutions for O(e^). The 
ODEs (|l^) are linear with inhomogeneous terms -Fj,ofc- Hence we construct homo- 
geneous solutions first. They are linear combinations of two independent solutions 
vi{t) and V2(t) written as 

vi{t) = -sinh(t)sech^(t), (19) 

3 13 
V2{t) = -sech(t) - -cosh(t) - -tsinh(t)sech2(t). (20) 

Note that vi(t) and V2(t) are odd and even functions of t, respectively and Wron- 
skian W = 1. One can construct solutions to the equations (|1^) under the 
conditions (|T5D and ([T8| ) by using the equations (|l^) and (pOD ; 



g,,oi(t) = isech3(t)-^^^^sech(t) + ^sinh(t)sech2(t). (21) 

In this way one can construct even solutions qj^iit), I > 2 successively. 

This result leads to no homoclinic bifurcation, which is inconsistent with the 
fact that the symplectic map (^) becomes chaotic when e 7^ 0. This indicates 
a breakdown of the apphcation of the regular perturbation method (|12]). The 
breakdown can be clearly seen near the singular points of the outer solutions (^) 
and (ETl). We will discuss the breakdown in the next section. 



3 The breakdown of the outer expansions 



In this section we analyze the singular points of the outer solutions (|TB|) and (PT]), 
extending the domain to complex t-plane. 

The unperturbed solutions (|16]) have first order poles at t = 7rz/2+7rm, -n, G Z. 
And the solutions for O(e^) ( [2T| ) have third order poles at the same points. One 
can show that the solutions for 0(e^') have the {21 + l)-th order poles (in general. 
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with ramification) at t = 7ci/2 + vrm, n G Z, i.e., for example, magnitude of 
qjfii{t) in tlie neighborhood of t = 7ri/2 are given by 



{t-ni/2) 



21+1 



o 



TT 

t i 

2 



(22) 



where aj^ are pure imaginary constants, for example, ai^ = 02,0 = —i, = 
^2,1 = V3- By substituting the expressions (j2^) into the expansions (|12|) one 
obtain 



21 



1=0 

aj,ie 



V- 



T 1 + 



TT 

'--2' 



(23) 



From the equations ( psD one can see that all the terms in the expansions (0) 
give contributions of the same magnitude in |t — 7ri/2| ~ e and that the regular 
perturbation method breaks down there. In the next section, we analyze the 
region near the singular points by blowing up there. The region are called the 
inner region, whereas the region in which the regular perturbation method is 
valid (t — > — 00) is called the outer region. 



4 The inner equations 

In the outer region we have obtained the even solution of the sub-manifold. This 
solution is valid in asymptotic of t — —00, but near the singular points the 
asymptotic expansions (|T2p break down and one can't construct analytic contin- 
uation of qj,u{t, e) from the left half of t-plane to the right one. To construct the 
analytic continuation and capture odd parts of the manifold, we shall blow up 
near one of the singular points, t = 7ii/2 as 

t-'^i = ez. (24) 

In addition, we transform qj{t,e) to ^j{z,e) in order that the equations to solve 
possess no negative powers of e as 

eqj{t,e) = <l>j{z,e). (25) 

By substituting the equations ( p^ and (^) into the equations (0), we obtain 

A2$,(z,e) = -2$,(z,e)=^ + e'<l>,(2,e) -2/^e'$,(2,e)<l>fe(^,e)^ (26) 

where A^$j(z, e) = ^j{z + l,e) - 2^j{z,e) + ^j{z - l,e). The second order 
difference equations (B^) are called the inner equations. 
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In this paper we try to construct the leading term of ^j{z, e), hence we restrict 
our attention only to the leading term of the inner equations, i.e., 

A^^z) = -2$(z)3. (27) 

We also call the equation (|27D the inner equation^. 



The formal solution of the inner equation and 
Borel transformation 



One can find the formal solution to the equation (l27| ) around z = cxd in the power 
series of 1/z as 

= (28) 



1=0 



where a/ =^ ai^i = 02,/, / = 0, 1, ■ ■ ■ (see the equations (^)). Because <l'o(-2) is a 
solution to the equation (pT]) with boundary condition ^o{z) ^ for \z\ —>■ 00, 



the coefficients a^, / = 0, 1, ■ ■ ■ satisfy the following estimation|[16||: 

ai ~ (-l)'+i(2/ + 3)!(27r)-2'-2 (/^oo). (29) 

Although the radius of convergence of the right hand side of the equation (^) 
is obviously 0, ^o{z) is an asymptotic expansion of ^u{z) for Re 2; < and of 
^s{z) for Ke z > 0. Additionary ^o{z) is Borel summable, i.e., the formal Borel 
transformation (the inverse Laplace transformation) of ^o{z) 

V{p) ''^ L-'%{^) (30) 
00 

has a finite radius (= 2it) of convergence and the solution $(-2) = LV{p) is known 
as the Borel sum of ^o{z). 

We hence apply the formal Borel transformation to the inner equation (p7|) 
and obtain 

2(cosh(p) - l)V(p) = -2 [V *V *V ]{p), (31) 

where V{p) has been defined as the equation (|30D and convolution [ V"*iy *F](p) 
is defined as 

[V*W*Y]{p) = [[V *W ]*Y ]{p) 



[V*[W*Y] ]{p) 

rv{T) r " w{p-T - x)Y{x)dxdT. 

Jo Jo 



^In the leading equation of the inner equation is caUed the truncated inner equation. 
Here we cah the equation ( p7[ ) the inner equation for simphcity. 



The solution to the integral equation ( pT] ) possesses poles only at p = 27iik, k G 
Z\{O}0. Note that these singular points are determined by the kinetic terms of 
the generating function, hence in many physical problems this type of singular 
points appears. 

One can define the unstable sub- manifold ^u{z) and the stable one ^s{z) by 
using the Laplace transformation as 

$,(2) = / e^PW{p)dp , |argz| < tt, (32) 
Jo 

r — OD 

$„(z) = / e~PW{p)dp , -271 < argz < 0, (33) 
Jo 

respectively. Now let the difference between ^s{z) and ^^(-z) be ^±{z) when 
±Im z > 0, i.e., 

$±(z) = ^s{z) - ^uiz) (±Im z>0). (34) 
From the equations (|3^) and (^) one obtains 

zIm<l>„(-zC) = -i$-(-<), (35) 

where ( is an arbitrary positive real number (see Appendix 0). 
Furthermore from the equations (P^), ( |5B| ) and one obtains 

$_(z) = f e-f"V{p)dp {-7i< aigz <0). (36) 

The integral path 7 is shown in figure 0. $-(-2) is given as the sum of the residues 
of the poles on the positive imaginary semi-axis, i.e., 

00 

$_(2) = 2niJ2^es[e-P'V{p)]p = 2mk]. (37) 
fc=i 

Because of |e~^'^'^| ^ |e~^'^*^| ^ ■ ■ • in Im 2 < 0, we shall only take account 
of the first term of the right hand side of the equation (|37|), i.e.. 



$_(2) ~ 27iiRes[e-P^V{p);p = 2m]. 

From the equation (^) one can realize that $u(^) and ^s{z) possess the term 
in proportion to e"^'^*'^^, hence we look for solutions to the equation (pTl) with 
the form 

$(z) = 5:<l'K^)e-2"*'^ (3g) 
/=o 
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Figure 1: The integral path 7 on complex p plane is shown. One can take 7 as 
enclosing the positive imaginary semi-axis. 



By substituting the equation (|38D into the equation (|27D , equations are succes- 
sively obtained for each power of e"^'^*^; 



0(e°) 



A2|.o(z) 



-2<^>o(^)^ 



(39) 
(40) 



The formal solution to the equation (^) has been obtained as the equation 
(pHI). In order to find the leading term of a solution to the equation (^), we 
replace $o(-2) by its leading term —i/z. This yields 

A'^,{z) = ^^,{z). (41) 



The equation (|l|) possesses two formal solutions •I'f 



z), $° (z) whose leading 
^ respectively. Because the odd solution ^'f'^{z) dominates the 
even solution $^'"^"(2;), we take account of only ^°'^'^{z), i.e., we replace ^i{z) as 



terms are z ^, z , 



odd/ 
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One can continue the process and obtain the formal solutions of / > 2. In 
this paper, however, we restrict our attention only to the first two terms of the 
equation (p^), $o(-2) and $1(2). 

The odd solution <l>f^( z) possesses third order poles at z = 00. Hence it 
is easily proved that V{p) possesses fourth order poles at p = 2'7iik, /c G Z by 
properties of Laplace transformation, i.e.. 



V{p) 

From the equation 



dk 

{p - 2mky 



+ 



dk : const, k E Z. (42) 



we can obtain ^~{z) as 

1 



2iTi lim 



2^27ri 3! (ip3 



[{p - 27riye-P'V{p)] 



where c ( = —27Tidi/3\ ) is a Stokes constant to be evaluated. There have been 
no analytical methods to evaluate the Stokes constant. However it is possible to 
estimate numerically by using the same algorithm as used in [p|-p!2[]. The result 
of our calculation shows c ~ —796. 

By using the equation ( P5| ) one can let the coefficient of bz^e~'^'"^^, which is 
the leading term of $^^1(2;), as 



TT 

, — TT < arg2; < 

C TT ^ 

- , arg2; = 

2 ' Z 2 
-c , — — < arg2; < 



(43) 



approximately. 



6 Matching of the solutions 

In this section we match the inner solutions to the outer solutions and construct 
the unstable sub-manifold which is valid on the right half of t-plane. 

First we give the exponential expansions of the outer variables. Following the 
inner expansion (|38|) , the refined outer expansions are given by 



Qj (t, e) = gj- o(t, e) + qj,i{t, e)e + g^, 2(t, e)e + ■ ■ ■ . (44) 



Note that the exponential expansion ( ^4|) is singular with respect to e, i.e., e 
has essential singularity at e = 0. 
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The leading term of the expansion (^) has been obtained in section |^ as 
qj,o{t, e) = gj,oo(t) + e^gj,oi(i) + O(e^), 

where 

qj,oo(t) = sech(t), (45) 
gj-oi(t) = ^ sech^(t) - ^ ^^^^ sech(t) + sinh(t)sech^(t). (46) 

Next we construct the first exponential terms gj,i(t,e), which correspond to 
0(e~^*). Substituting the expansions ( ^4]) into the equations (^ gives 

A^gj,i(t, e) = (1 - 6gj,o(^, e)^)gj,i(t, e) - 2e^Kqj^o{t, e)qkfl{t, ef . 

To estimate the leading order of qj^wit), which is the leading term of gj,i(t, e), 
we replace qj^it^e) with their leading term gj^oo(^)- By these replacements, we 
obtain the differential equations of qj^ioit) as 

= (l-6g,oo(t)^)g,,io(t). (47) 



The general solutions of the equations (|47|) are given as 

qj,10{t) = Cj-ifi(t) + Cj-2f2(t), 

where Cj^i and Cj^2 are arbitrary constants, and Vi{t) and V2{t) are given by the 
equations ([T9|) and (^Ol). 

We rewrite functions fi(t) and V2{t) in the neighborhood of the singular point 
t = 7Ti/2. If we introduce 6 = t — TTi/2, the equations (|19[) and (pO|) imply 

, , cosh (5) , , 

3 

--i ^[3sinh(5) + smh\S) - 3Scosh{S)]. (49) 

2 sinh (o) 

One can see that the odd function Vi{t) becomes an even function of 6 and the 
even function V2{t) has both odd and even parts. The leading term of vi{t) around 
5 = is of the order 0(5"^) = 0(e^^z~^). It corresponds to the solution ^'^"'^'^{z) 
of equation (0), which is sub-dominant to the other solution $°'^'^(z). The odd 
part of V2{t), which is in the square bracket, has the Taylor expansion as 

3sinh(5) + sinh3(5) - 35cosh(5) 

= 0-6 + 0-6'^ + -6^ + 0(6^). 
5 
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The coefficient of tlie square bracket is of the order 0(5^^). These show that the 
odd part of the equation (^9[) is of the order 0{6^) (= 0{e^z^)), corresponding to 
|,odd^^^_ We shall therefore match qj,ioit) to ^i{z) = ^^''■{z), i.e., qjAoit) is given 
by the odd part of V2(t) as 



V2{t) - -TTWiit) 



(50) 



5 

= -lc,-22eV(l + 0(52)). 



Furthermore the discussion in section ^ shows 

i Im ~ hz^e-"^^" 

where h has been determined as the equation p3|). The matching qj,io(t) and 
^odd^^^ for t approaching 7ri/2 with arg2; = — 7r/2 is given below: 



egj-io(t)exp 



'-cz^e-^^'\ 



hence, 



C7,2 



5 . c 

--«-exp I — - 



TT 

e 



(51) 



where c G R, therefore, c,-^2 ^ ^R- Another constant Cj^i is given as 



15 c 
— TT— exp 
8 ^ 



TT 

e 



(see equation (|50|) and (^Tj)). 

Finally, the leading term of the imaginary part of the sub-manifold qj,io{t) is 
given as 

gi,lo(^) = Cj^iVi{t) + Cj^2V2{t) 

1 1 o 

= --^q,2 . , o/rx [3sinh((5) + sinh (5) - 35cosh((5)] 
/ smh (o) 



9^ 
2 



3(t — — z)sinh(t)sech [t) — 3sech(t) + cosh(t) 
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7 Results 



We have constructed the solutions which are vahd in the neighborhood of the 
singular point t = 111/2. The same arguments, applied to the complex conjugated 
singular point t = —ni/2, yield the complex conjugated solutions. We take the 
two singular points, which are nearest to the real axis, into consideration and 
restore the obtained solutions to real functions. Taking the average of these two 
solutions, we obtain 



Qj (t, e) = qj^o{t, e) + Re 



^gj-z(t,e)exp 
i=i 



27ri/ 



-t 



In this paper we approximate the sub-manifold by taking the term / = 1. This 
term is written as 



Re 



gj-i(t,e)exp 



2mt 



Re 



gj-io(t)exp 



2mt\ 



where 



M{t,e) 



Cifi(t)cOS lC2V2it)sm 

e e 
M{t,e)B{t,e) 



3 , , 2nt , , 27ct 
-7rfi(t)cos i;2(^)sin . 



(52) 
(53) 



The switching function S{t,e), which stands for Stokes phenomenon, is approxi- 
mately given as 



S{t,e) 



(t < 0) 

1 (t = 0) 

2 (t > 0) 



(54) 



(see the equations (^31) and (|5^)). 

By using these functions, the parametrization of the unstable sub-manifold is 



g,-,(t, e) = g,-oo(t) + e\j,oi{t) + M{t, e)B{t, e), 
and the momentum variables Pj,u{t, e) is given as 



(55) 



(56) 



(see the equation (|^)). On the other hand, the parametrization of the stable 
sub-manifold is given as qj,s(t, e) = qj^u{—t, e) and Pj,s(t, e) = —pj,u{—t + e, e). 
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By using the equations ( |56D one can observe the manifold on the phase space. 
However it is impossible to display a 4-dimensional phase space graphically, hence 
we observe its projections to {qi,pi) and (^2,^2) plane. From equations (|55D , (^) 
and ( ^61) the projections to and (^2,^2) plane are same. Figure ||,(a) shows 

that global shapes of the stable and unstable sub-manifolds. In figure ^(b)-(d) we 
magnify the neighborhood of the fixed point. One can clearly see the homoclinic 
structure from figure |],(d). 

Next we compare the approximate solution (p5D and (|56|) with numerical 
trajectories. 20 numerical trajectories are plotted on figure ^. This figure shows 




Figure 3: The comparison between the approximate solution of the unstable sub- 
manifold and 20 numerical trajectories. The parameters are same as figure ^ 
The initial values of the numerical trajectories are given on the solution between 
t = 14.0 and t = 14.2 with equal intervals. 



that the approximate solution agree well with the numerical trajectories. The 
initial values of the numerical trajectories are given on the approximate solutions 
between t = 14.0 and t = 14.2 with equal intervals, i.e., t = 14.0, t = 14.005, 
■ ■ ■, t = 14.195. In other values of e or k also one can show that the approximate 
solution agree well with numerical solutions. 



8 Summary and discussions 

We extended the method "asymptotic expansions beyond all orders" for the 4- 
dimensional symplectic mapping straightforwardly and successfully obtained 
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the explicit functional approximation of the 1- dimensional sub-manifolds of the 
stable and unstable manifolds (^). These asymptotic solutions well approximate 
numerical solutions in the neighborhood of the hyperbolic fixed point, which is 
important to understand phase space's structure of Hamiltonian systems. We 
therefore consider that by this extension one can analyze global structure of 
phase space of Hamiltonian systems with three degree of freedom. 

The oscillating parts M{t, e) and B{t, e) do not contain the coupling con- 
stant K (see the equations (|52|) and (|53|)). This indicates that the homoclinic 
structure is a direct product of 2-dimensional symplectic mappings. In other 
words, the model we analyzed shows only quantitative characteristic properties 
of high-dimensionality. This is based on the fact that the inner equation (^) 
is decoupled. From the order of the singular points of the outer solutions, we 
find the condition that the coupling term dominates the nonlinear terms. The 
condition is written as 

J «! + 0:2 < 7-1-4 ■•• decouple , . 

) ai + a2 > 7 -l- 4 ■ ■ ■ coupling term dominates ^ ' 

(see the map (|lD and appendix |B|) . When the second condition of the conditions 
( |57D is satisfied, the coupling term dominantly contributes the inner equations 
and another type of rescaling (see equations ( p4|) and (|25|) ) is required. Hence we 
consider that the models whose coupling terms satisfy the second condition make 
not only quantitative but also qualitative differences. However it is impossible to 
treat such models with the same way. We would like to improve this method in 
order to extend to such models in a future work. 

In the present work we have restricted our attention to the particular 1- 
dimensional sub-manifolds of the stable and unstable manifolds. And we have 
shown that the sub-manifolds oscillate in the vicinity of the origin. Because the 
unstable sub-manifold is contained in the unstable manifold, it is clear that the 
unstable manifold itself oscillates. We now try to construct explicit functional 



approximations of the whole 2-dimensional stable and unstable manifold |jT7|. 

In the present model there exists a 4-dimensional fully hyperbolic fixed point 
(a direct product of two 2-dimensional hyperbolic fixed points). In a future 
work we will extend this method to 4-dimensional symplectic mappings with a 
hyperbolic X elliptic type fixed point. 



Acknowledgements 

We would like to thank K. Nozaki for helpful discussions and suggestions. 



References 



16 



[1] Guckenheimer J and Holmes P 1983 Nonlinear oscillations, dynamical sys- 
tems, and bifurcations of vector fields (New York: Springer) 

[2] Neishtadt A I, 1984 J. Appl. Math. Mech. 48(2) 133 

[3] Lazutkin V F, Schachmannski I G and Tabanov M B 1989 Physica 40D 235 



[4: 

[5] 
[6] 
[7 
[8 
[9 
[10 

[11 
[12 

[13 

[14 

[15 

[16 

[17 



A The imaginary part of ^u{z] 

In this appendix we prove equation (p5|). From equations ( |32D and (|33| ) we obtain 
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K{-z) = / eP'Vip)dp 
Jo 

POO 

= / e-P'Vi-p){-dp) 
Jo 
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and 

^s{z) = -^s{z). 

Assume C £ R- One can easily prove that 

$„(-<) = -$,«) = 2>,(-<). 
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Furthermore assume C > 0, and and are real- valued on the negative imag- 
inary semi axis, one can obtain 



(5^ 



By using equation ( ^8|) and equation (|59[) , we can say that 

$,(2;) = <i>u{z) {z e iR). 

Hence Re = Re ^u{z) on the imaginary axis: i.e., 

e G iR). 



Particularly when C > 0, from equation ([58|) one can obtain 



= -2Um$,(-zC). 



Thus 



(59) 



iIm<l>„(-zC) = --$-(-<). 

B The order of the couphng term 

In this appendix we study the order of the coupling terms in the map (|l]). The 
coupling term is generated by e'^J = (-'^ qi%,^2%^ C(j ^ N, j = 1, 2. 

Put A = ai + a2. First let us suppose 7 = 2. The outer solutions of O(e^) 
are written as 



9i,oi(t) 



C{t-A) 



1 7 

-sech-*(t) sech(t) 

3 ^ ^ 24 ^ ' 

+ ^ sinh(t)sech^(t) + KajC{t; A), 

-sech^(t) 
3 ^ ' 

1 (A -4)!! 
I (A -5)!! 

xE l1~!~o5;; ^(^;--4+i+2fc) (a>4) 



(A = 3) 



(60) 



(61) 



^-^ {A-2-2k)\l 
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where (2ri)!! 
f-l)!! = l, ( 



-- 2n-(2n-2) 
-2)!! = (-3)!! 

F{t-n) = 
gd(t) = 



2, (2n-l)!! = (2n-l)-(2n-3)---l,n G N, 0!! 

— 1, and F{t; n),n E Z is defined as 



f cosh"(t) {n > 1) 

j sinh(t)sech2(t)gd(t) {n = 0) 

tan~^(sinh(i(;)) 



(62) 



and 



M 



def 



^1 

- - 1 
2j 

2 

A+1 



(63) 



• A : even 
■ A: odd 



From equations (|5Il|)-(p3|) one sees that when A > 6 the outer solutions 
Qj,oi{t), J = 1; 2 possess singularities of the order A— 3 (> 3) at t = 7ii/2+7iik, k G 
Z (in general, with ramification). On the other hand, when A < Q, the order of 
the singularities is 3. One can easily prove that the outer solutions of 0(e^'^) 
possess singularities of the order k{A — 4) + 1 or 2A; + 1 when A > Q or A < Q, 
respectively. Hence if A > 6, the magnitude of qj,ok{t) in the neighborhood of 
t = TTi/2 are given by 



Clj,k 



J'j,k 



{t-m/2f^^'^^' 
C, j = l,2, k>0. 



O 



t 



IT 



(64) 



From equations (0) one can realize that all the terms in the formal expansions 
(|I2[) give contributions of the same order in \t — ni/2\ e^,S = We must 

therefore blow up near the one of the singular points, t = 7rz/2 as 



TT 



(compare equation ( plf ) and the discussions in section y). Note that if A > 6, 
then 5 < 1, hence > e, that is, the inner region is wider with A > 6 than with 
A<Q. 

Now we show that when 7 = 2 and A > 6, the rescaling factor is different 
from that of A < 6. Hence we call the models with A > 6 or A < 6 strongly 
or weakly coupled model, respectively. In the same way, one can prove that the 
condition with which the rescaling factor is different from 2-dimensional mappings 
(storongly coupled model) is written as 



tti + ^2 > 7 + 4 



(65) 



for general 7 G N. 
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